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ABSTRACT 

The light curves of some luminous supernovae are suspected to be powered by the spindown energy 
of a rapidly rotating magnetar. Here we describe a possible signature of the central engine: a burst of 
shock breakout emission occurring several days after the supernova explosion. The energy input from 
the magnetar inflates a high-pressure bubble that drives a shock through the pre-exploded supernova 
ejecta. If the magnetar is powerful enough, that shock will near the ejecta surface and become 
radiative. At the time of shock breakout, the ejecta will have expanded to a large radius (~ 10®^ cm) 
so that the radiation released is at optical/ultraviolet wavelengths {Tgs « 20, 000 K) and lasts for 
several days. The luminosity and timescale of this magnetar driven shock breakout are similar to 
the first peak observed recently in the double-peaked light curve of SN-LSQ14BDQ. However, for a 
large region of model parameter space, the breakout emission is predicted to be dimmer than the 
diffusive luminosity from direct magnetar heating. A distinct double peaked light curve may therefore 
only be conspicuous if thermal heating from the magnetar is suppressed at early times. We describe 
how such a delay in heating may naturally result from inefficient dissipation and thermalization of 
the pulsar wind magnetic energy. Without such suppression, the breakout may only be noticeable 
as a small bump or kink in the early luminosity or color evolution, or as a small but abrupt rise 
in the photospheric velocity. A similar breakout signature may accompany other central engines in 
supernovae, such as a black hole accreting fallback material. 


1. INTRODUCTION 

Optical surveys are finding a growing number of bril¬ 
liant, though rare, explosive transients, some 10 — 100 


dQuimby et al.|2007 

Smith et al.||2007 

Ofek et al.||2007 

Barbary et al.|2U091 

Pastorello et al.|20] 


Cuimby et al. 

2011 

Chomiuk et a 

1120111 Cal-Yam||2012 

Howell et al. 

2013 

inserra et al.|20131|iNicholl et al.|2014 

IPapadopou- 

los et 

al. 2015). The mechanism generating these enor- 


power ordinary supernova light curves - the diffusion of 
shock deposited thermal energy, or heating by radioac¬ 
tive ^®Ni - appear incapable of reproducing the observed 
rise time and peak brightness of many of the super- 
luminous supernovae (SLSNe). 

Two classes of models are frequently invoked to explain 
the light curves of SLSNe. The first involves interaction 
of the supernova ejecta with an extended circumstellar 
medium (CSM). If interaction occurs at a location where 
the ejecta is translucent (radii ^ 10®® cm) the thermal- 


ized kinetic energy can be radiated efficiently ( Wooslgy 
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Moriya et al.||201ip. In tne second class of models, 


the ejecta from a seemingly ordinary supernova explo¬ 
sion is continuously reheated via energy injection from a 
long-lived central engine, either a rapidly rotating, highly 
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the class of Type H SLSNe, narrow (^ 10 — 100 km s“®) 
hydrogen Balmer lines are often seen in emission, a sig¬ 
nature of interaction with a slow moving CSM. In the 
Type I SLSNe, on the other hand, no spectroscopic indi¬ 
cations of interaction are apparent and the line features 
only indicate rapidly moving (^ 10,000—15,000 km s“®) 
material. The magnetar powered model of SLSNe has 
been successful in fitting the light curves, colors, photo- 
spheric velocity evoluti on, and gross spectral features of 
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uing for lOO’s of days after the explosion, which has been 
claimed to be indic ative of persistent magnetar heating 
( Inserra et al.|[20l3 ). 

Additional observational tests are needed to validate 
and discriminate models of SLSNe. Such an opportu¬ 
nity may have arise n with the well-sam pled photometry 
of SN-LSQ14BDQ dNicholl et al.||2015D, a Type I SLSN 
with a double-peaked light curve. T he luminosity of SN- 
LSQ14BDQ rose to an early maximum in « 5 days; then, 
after a brief decline, the light curve rose again to an even 
brighter peak (« 2 x 10®® ergs s“®) by 50 days later. A 
similar double-peaked morphology had al ready been seen 
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tral engine heats the supernova ejecta in two physically 
and spatially distinct ways. A thermalized pulsar wind 
heats the ejecta directly at its base, powering a luminos¬ 
ity that diffuses out on timescales of weeks or months. 
At the same time, the pulsar wind dynamically affects 
the ejecta, inflating a high-pressure bubble that drives 
shock heating at larger radii. If the magnetar is power¬ 
ful enough, that shock will near the ejecta surface and 
become radiative, producing an early burst of emission. 
The situation resembles, in some ways, shock break¬ 


out from a stellar explos ion (Klein & Chevalier 1978 


Matzner & McKee 19991, with a few key differences. 


h'irst, a magnetar driven shock propagates through a 
moving medium; the shock will be weaker, and when 
it does emerge, the ejecta surface will have expanded 
by several orders of magnitude in radius. The resulting 
emission will last longer and be at longer wavelengths 
(optical/UV) than the brief x-ray burst that accompanies 
ordinary supernova shock breakout. Second, the shock 
does not necessarily die once it becomes radiative; as 
long as the magnetar continues to inject energy, the shell 
can be driven faster than free-expansion and may release 
energy at its outer edge. 

The integrated heating from a magnetar driven shock 
amounts to only a few percent of the total pulsar wind 
energetics. However, shock heating occurs exterior to 
the bulk of the ejecta, and so can be radiated < 1 week 
after explosion, before most of the centrally thermalized 
energy has had time to diffuse out. Under certain cir¬ 
cumstances, the two heating mechanisms may produce 
two distinct emission maxima. We develop below an an¬ 
alytic description of “magnetar driven shock breakout”, 
and present toy light curve calculations that suggest that 
this effect provides an appealing explanation for double- 
peaked supernova light curves, and a means to constrain 
the magnetar model of SLSNe. 


2. DYNAMICS OF MAGNETAR DRIVEN SHOCKS 

Consider a pulsar with spin period P and magnetic 
field B. The total spin energy is 

An, « 2 X 1052p-2 ergs, (1) 

where P^s = P/^ ms and we adopt a neutron star mo¬ 
ment of inertia of J = 10'^® g cm^. The rate at which 
energy is input from the pulsar is, in the case of vacuum 
magnetic dipole spindown 




Fig. 1.— ID hydrodynamical calculation of the dynamical effect 
of a magnetar on the supernova ejecta. The ejecta had Msn = 
5 M 0 , and initially Esn = 10^^ ergs and a broken power law density 
profile (Eq. [^. Energy injection from a magnetar with Em = 
5 X 10^^ ergs and tm = 5 days was input according to Eq. in 
the inner few zones of the ejecta. The snapshot shown is at a time 
t = 15 days. The blue shading highlights the low density, high 
temperature magnetar bubble. The red shading shows the region 
heated by the magnetar driven shock. Radiation diffusion has not 
been included in this calculation. 


sity profile, with a shallow profi le in the inner region an d 
a steep one in the outer regions (Chevalier & Soker 1989). 
The transition occurs at a velocity coordinate 


vt = 


(4) 


In the inner ejecta (v < Vt) the density at a position r 
and time t is 


P{r, t) = Cp 




(5) 


The density profile in the outer ejecta {v > Vt) has the 
same form but a different exponent, p oc r”". The coef¬ 
ficients are given by 


(n-3)(3-,5) 
47 r(n — (5) 


( 6 ) 


Lm — 


Ernttu 




where the spindown timescale is 




c; n-2 p2 

on>l4 


days. 


( 2 ) 


(3) 


where P 14 = P/1 0^^ gauss and we assume, as in Kasen 
& Bildsten (2010), that the angle between the rotation 
axis and magnetic dipole is a = 45°. If the spindown en¬ 
ergy is thermaliz ed in the ejecta, it may power a super¬ 
nova light curve dBodenheimer fc Ostriker 1974 Gaffet 


1977 |Maeda et al.|2007| ' Kasen fc Bildsten|26lO[|Woosley 


2010 1 . The dynaniical effect of the energy injection, how- 


ever, is independent of whether it thermalizes or not; in 
either case the energy behaves as a 7 = 4/3 gas. 

We will assume that the supernova ejecta of mass Mgn 
is spherically symmetric and has a broken power law den- 


G = 


2(5-5)(n-5) 

(n-3)(3-(5) 


- 11/2 


(7) 


Ty pical values for core collap se supernovae are S = l,n = 
10 ( Chevalier fc Soker|1989 ), which we adopt as fiducial. 

For Am A'sn, the magnetar wind will significantly 
restructure the supernova ejecta, as illustrated in the 
hydrodynamical calculation of Figure The high pres¬ 
sure from central energy injection evacuates a cavity and 
sweeps ejecta into a thin shell. The shell moves faster 
than the local ejecta expansion velocity, and a radia¬ 
tion dominated shock of relative velocity « 2500 km s“^ 
forms. In a multi-dimensional calculation, Rayleigh- 
Taylor instabilities would b reak apart the shell a nd smear 
out the density peak (e.g., Blondin et al.j 2001), but the 
global structure would be qualitatively similar. In addi¬ 
tion, at late times, radiation diffusion from the inner hot 
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bubble will smear out the temperature discontinuity at 
the shock front. 


2.1. Time of Shock Emergence 

The magnetar driven shock will become radiative when 
the diffusion time to the ejecta surface, tdiff ~ ri?/c, 
becomes comparable to the elapsed time, t ^ R/v, or 
when the shock reaches an optical depth t cjv. In the 
outer ejecta, the optical depth from the surface inward 
to velocity coordinate v is 


where 
td = 


t{v) 



c ti 


1) Vt f^ 





-I 1/2 


VtC 


32 days. 


( 8 ) 

(9) 


is the effective di ffusion time i n a homologous ex- 


^sn,5 — 


panding medium (Arnett 1982). Here Mg, 
Msn/5MQ,Esn,5i = E^n/W^ crgs, and kq.i = 
k/O.I cm^ g“^ is the scaled opacity. The shock then 
becomes radiative (t = c/v) when it reaches a radius 


^bo = Vtt 


1 


td 


_^/n — 1 t 


2/(n-2) 


( 10 ) 


Eq. 10 assumes that r^o/t lies above the transition ve¬ 


locity coordinate Vt, which is true for t < tdj'ijn — 1 
10 days. 

To determine the time when the shock reaches the 
breakout radius rbo, we make the assumptio n that mass 
is swept up into a geometrically thin shell ([Os tri^ej^^Z 


GunnHioTT Chevalier|1984 Chevalier fc Fransson|199’2 r 
i'he momentum and energy equations describing the 
shell dynamics are 


M 


dvs 

dt 


= Aiirllph - p{v^ - 


( 11 ) 


d(47rr3pb) . 2d?’s , , 


(i)-L„(f), (12) 


where Mg, Ts, are the mass, radius and velocity of the 
shell, Vej = Ts/t is the ejecta velocity at radius rg, p is 
the preshock ejecta density ahead of the shell, and pb is 
the pressure in the magnetar inflated bubble. The Lgn 
term represents the rate at which thermalized magnetar 
wind energy radiatively diffuses out of the bubble. 

Assuming that the magnetar injects energy at a nearly 
constant rate, Lm = E^ftm and that diffusion losses can 
be ignored (Lgn ~ 0, appropriate for t <C td), the dynam¬ 
ical equ ations have self-similar powe r law solutions for rg 
and pb ( Chevalier fc Fransson||1992 | 

rs{t) = (13) 

where a = (6 — i5)/(5 — (5) and 

hr = Ctr (14) 


is the time it takes to the shell to propagate through the 
inner ejecta and reach the transition velocity Vt- The 



Fig. 2.— The different regimes of magnetar driven shock emer¬ 
gence, illustrating whether the shock becomes radiative in the in¬ 
ner ejecta (ubo < or outer ejecta (ubo > '^t) and whether the 
magnetar is “on” (t^o < ^m) or “off ” (tbo > ^m) at the time of 
emergence. The solid lines are determined from Eqs. |18| and |19| 
with a = 5/4 and n = 10. Breakout emission can occur in all re¬ 
gions, but will be brightest in region 2, when the shock is strongest. 


coefficient is 


Ctr — 


2(n-5)(9-2(5)(ll-25) 


(5-5)2(n-<5)(3-(5) 


(15) 


For u = 10,(5 = 1, we have a = 5/4 and Ctr = 2.2. The 
expression for itr has been previously derived in|Chevalier| 

(po^. 

'i’his self-similar power law solution only holds for 
times t < ttr when the shell remains in the inner ejecta. 
At later times, the shell moves into the steep outer ejecta 
and accelerates. The shock front will begin to move 
ahead of the shell, but the thin shell approximation ap¬ 
pears to remain reasonably valid (see Figure [^. The 
limiting behavior can be understood by noting that, at 
large radius, nearly all of the mass is swept into the shell. 
Mg —>■ Msn- Civen the low densities in the outer ejecta, 
the ram pressure term p(vs — Uej)^ can be neglected rela¬ 
tive to the bubble pressure, and the d yna mical equations 
give asymptotically a power law (Eq. 13) with a = 3/2. 

Although the break in the density profile complicates 
the dynamics, the shell radius can be reasonably approx¬ 
imated by rg oc where the exponent is in the range 
a « 1.2 — 1.5. If we assume that the shock becomes ra¬ 
diative in (or near) the outer ejecta, the time of shock 
breakout is found by setting rg equal to the breakout 
radius, rbo (Eq. 10), giving 


tbo 


.1-/3 

Str 


(n- l)/3/2 




(16) 


where 



2 

l){n-2) + 2' 


(17) 


The shock breakout time, tbo, is a weighted geometric 
mean of td and tm, where the appropriate value of jd 
depends on whether the shell has spent most of its time 
in the inner or outer regions of ejecta, i.e., whether tbo is 
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Fig. 3.— Numerical calculation of the shock breakout time, tbo? 
for various values of Em and tm, and standard ejecta parameters: 
Man = 5 Mq, Esn = 10 ^^ ergs, S = l,n = 10, K = 0.1 g cm^, 
which imply a diffusion time = 32 days. The magnetar energy 
injection rate was taken from Eq. [21 The solid line shows the 
estimated = im contour from Eq. 119| the dashed line shows the 
estimated = '^t contour from Eq. |18| 


much greater than or less than ttr- For tbo ttr, the self¬ 
similar value a « 5/4 applies, and /3 = 1/2 (for n = 10). 
For tbo 3> ttr, the asympotic value a = 3/2 is more 
accurate and /3 = 1/3. The distinguishing condition is 

tbo > ttr if (18) 

td 


To further complicate the dynamics, at times t > 
the energy input from the magnetar “shuts off’, i.e., 
drops towards zero. Following this, the bubble pressure 
decreases nearly adiabatically and the shell radius ap¬ 
proaches free expansion, Ts oc t. A shock breakout can 
still occur, but the emission will be less luminous because 
the shock speed declines for t> The brightest shock 
breakouts then occur when the shock becomes radiative 
before the magnetar shuts off, which requires a magnetar 
energy 


tbo < tm if E^ > O.TAs; 


td\ 

tmj 


/3/(l-/3) 


(19) 


The two conditions (Eqs.[l^and 19) divide the E^-t^ 
parameter space of shock emergence into four regimes, 
illustrated in Figure The partitioning is only sugges¬ 
tive, as real ejecta density profiles are more complex than 
a broken power law, and the magnetar energy deposition 
does not shut off sharply, but rather follows a smooth 
function of the form Eq. The most luminous break¬ 
outs occur when the shock is being continuously driven 
through the steep outer ejecta (region 2), which happens 
when the magnetar is energetic and/or long lasting. On 
the other extreme, for Em ^ Egn or ^ the shock 
will stall out before being revealed and no prominent 
shock breakout signature is expected.. 

To better determine the time of shock emergence, we 
num eric ally integrated the dynamical equations (Eqs. 11 
and [T^ to follow the shell location and find when it 
reached the t = cjv surface. This calculation did not 
assume Lgn = 0, but rather included approximate ra¬ 
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Fig. 4. — Numerical calculation of the shock strength parameter 
rj = (ds — Ve^)/vQ^ at the time of shock breakout for various values 
of Em and tm, and standard ejecta parameters: Mgn = 5 Mq, 
Esn = 10^^ ergs, 5 = l,n = 10, k = 0.1 cm^ g~^- The left 
panel assumes a constant magnetar energy injection rate Lm = 
Em/tm for t < tm and Lm = 0 afterwards. The right panel uses 
the continuous injection rate of Eq. The solid line shows the 
estimated fbo = contour from Eq. |W| the dashed line shows the 
estimated ribo = contour from Eq. |18| 

diative loses via the method described in the Appendix. 
Figure shows this calculation of tbo for various values 
of Em and tm, and standard ejecta parameters. When we 
carry out the numerical integration using Em constant for 
t < tm, and zero afterwards, the resulting tbo agrees well 
with the analytic expression Eq. For the more real¬ 
istic case of continuous energy injection (given by Eq. ^ 
the analytic result underestimates tbo by about 20%. In 
this case, a better estimate is achieved by multiplying tm 
in Eq. [^by a factor « 1.5 to account for the non-zero 
magnetar energy injection at later times. 


2.2. Shock Heating 

We next estimate the local heating from the magnetar 
driven shock, which will set the luminosity of breakout 
when the shock emerges. The rate at which energy is 
dissipated at the shock is 

esh(t) = 47rr2^(u, - Uej)^ = (20) 

where in the second expression we have parameterized 
the shock velocity as 


g{t) = 


Vs - Vej 


rs\dt t ) 


( 21 ) 


The shock strength parameter rj will play an important 
role in determining the radiated luminosity. If the shell 
radius obeys a powerlaw, Ts oc then rj = a — 1. In 
the inner ejecta, a approaches the self-similar value a = 
(6 — (5)/(5 — (5) and r] = 1/(5 — S), or t] = 1/4 for (5 = 1. 
In other words, the shell moves 25% faster than the local 
ejecta velocity. For t 3> Er, when the shock propagates 
deep into the steep outer ejecta and accelerates, a —)■ 3/2 
and T] —1/2. 

Using our numerical integration of the dynamical equa¬ 
tions discussed in Section we calculated the value of n 
at the time of shock emergence. The left side of Figure® 
shows the numerical determination of rj for the simple 
case where Em = E^/tm is constant for t < tm, then 
immediately drops to Em = 0. The behavior follows an¬ 
alytical expectations: in region 1 (shock emergence in the 
inner ejecta, magnetar on) 77 = 0.25; in region 2 (shock 
emergence in the outer ejecta, magnetar on) rj > 0.25 and 
increases with increasing magnetar energy, approaching 
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Fig. 5. — Calculation of the local shock heating rate (Eq. |20[ | over 
time for standard ejecta parameters Msn = 5 Mq, Esn = 10^ ergs, 
5 = Fn = 10, and continuous magnetar energy injection given by 
Eq.l^ The legend gives the value of the magnetar energy (in units 
of 10^^ ergs) and the spindown time (in days). A power-law break 
in the heating rate is seen when the shock enters the steep density 
profile of the outer ejecta. The luminosity of shock breakout is 
approximately the local heating rate at the time of breakout. 


a maximum value rj = 0.5. For region 3 (shock emer¬ 
gence in the outer ejecta, magnetar off) rj declines as tm 
decreases, illustrating the progressive weakening of the 
shock following magnetar shut off. 

The right panel of Figure shows the behavior of rj in 
a calculation using continuous magnetar energy injection 
given by Eq. Similar trends with Em and tm are seen, 
but the values of rj are generally lower, as the magnetar 
energy deposition is spread out over a longer timescale, 
with Lm always less than Em/tm- 

Figure shows the time dependence of the local shock 
heating rate (Eq.[20|) for continuous magnetar energy in¬ 
jection. The luminosity of the shock breakout pulse will 
approximately equal the local heating rate at the time of 
breakout. The local heating falls off with time roughly 
like a power-law, due to the progressive decrease in the 
preshock density. A break in the local heating rate is 
seen when the shock enters the steep density profile of 
the outer layers {v > Ut). The shock is weakened at times 
t > tm by the decrease in magnetar energy injection, and 
further at times t > hy radiation diffusion through 
the ejecta shell, which depressurizes the magnetar bub¬ 
ble. This latter effect is included in the calculation by 


including a non-zero diffusion term Egn in Eq. 12 


3. OBSERVATIONAL CONSEQUENCES 
3.1. Properties of Shock Breakout 

To estimate the peak luminosity resulting from shock 
heating, we can use “Arnett’s law” (Arnett 19821 which 


states that, for any general heating source, the luminosity 
at the light curve peak is equal to the instantaneous rate 
of energy deposition at that time. We thus determi ne t he 
peak luminosity by evaluating the heating rate Eq. 20 at 
the shock breakout time, tbo 


T MsnVf 3 f Vho 

Ebo = 27 rCp—- r] - 

tho \ Vt 


— n-|-5 


( 22 ) 


Shock breakout occurs at a velocity coordinate Vho = 


Tho/tho determined from Eqs. 10 and 16 


Vho « l.li^t 


Am 
i~ E 




(23) 


Evaluating Eq. at the velocity coordinate Ubo then 
gives the peak luminosity 


Ebo « 


. Aa 
td 


td A„ 


(„-8)(l-/3) 


(24) 


For standard parameters (/3 = 1/2, n = 10) the exponent 
of the term in brackets is only 1 / 8 , and we see that the 
peak luminosity depends most sensitively on the shock 
strength parameter, 77 . The luminosity also depends on 
Asn, as this sets the ejecta expansion velocity and hence 
the size of the remnant at the time of breakout. 

The spectrum of the breakout emission can be approx¬ 
imated by a quasi-blackbody with an effective tempera¬ 
ture, Teff, determined by Ebo = ^TcrpCrT^g, where r-p is 
the photospheric radius defined by the t = 1 surface. 
From Eq. 

rp«1.27;ttd(ibo/td)(”-"^/("-'). (25) 

Plugging fiducial parameters (n = 10, <5 = 1,/3 = 1/2) 
into Eqs. 24 and 25 we arrive at analytic estimates 
of the time, luminosity, photospheric radius, and effec¬ 
tive temperature at the time of a magnetar driven shock 
breakout 


^bo ' 

Ash ’ 
rp! 
Aeff ' 


(26) 


A 0.8 days 

.2.1 X 10^3 ergs st^7) 

.7.6 X cm (28) 

.15,500 


where As„, 5 i = Am/10®^ ergs, fm .5 = im/5 days, 770.2 = 
77 / 0 . 2 . In most stellar shock breakout events, the dura¬ 
tion of the burst is set by the light crossing time of the 
ejecta, however in the present case the ejecta is extended 
and the timescale is set by the diffusio n time, tbo (see 


e.g., Chevalier h Irwin]|2011 Piro||2013 ). 


The analytic results are only approximate; as an 
improved estimate of the peak luminosity, we used our 
numerical integration of the shell evolution to evaluate 
the heating rate at the time of shock emergence. 
Figure shows the numerical results in the E^-t^ 
parameter space. The luminosity increases with E^ as 
higher magnetar energy drives a stronger shock (greater 
77 ). For low values of tm, the the luminosity drops due 
to the decline of 77 following magnetar shut off. 


3.2. Approximate Light Curves 

The light curve of a magnetar powered supernova will 
be the sum of the emission from shock breakout and the 
diffusive luminosity from central magnetar heating. A 
first peak will only be distinguishable when the shock 
breakout luminosity is comparable to or greater than the 
diffusive luminosity at that time. To model t he com- 


posite light curve, we used a one-zone formalism (Arnett 
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Fig. 6.— Same as Figure[^but showing the calculated luminosity 
at the time of shock breakout. 


19821 to calculate approximate emission from each mech¬ 


anism, then added the results. The method is described 
in the Appendix. 

The left panel of Figure shows a composite model 
bolometric light curve. For reasonable magnetar 
and ejecta parameters {E^ = 5 x 10®^ ergs,ti„ = 
2 days, Mgn = 5 Mq, Agn = 10^^ ergs), the shock 
breakout emission is not dominant and produces only 
a kink in the early light curve rise. The lack of a promi¬ 
nent first peak is consistent with the grey radiation - 
hydrodynamical calculations of|Kasen & Bildsten (2010). 
Despite the lack of a distinct light curve bump, the break- 
out may still be detectable by the shift in brightness 
and colors at the time of breakout, or by a sudden, 
but small (~ 1000 km s“^) increase in the line veloci¬ 
ties when the photosphere recedes through the region of 
non-monotonic velocity. 

To see a clear double-peaked light curve requires ei¬ 
ther a very bright shock breakout, or a slowly evolving 
diffusive light curve, properties that are only realized in 
certain regions of parameter space. Figure]^ shows that 
increasing the ejecta mass delays the diffusive light curve, 
making the shock breakout peak more prominent. In¬ 
creasing the kinetic energy of the supernova explosion 
leads to a larger remnant and brighter shock emission, 
which also clarifies the double-peaked structure. In ad¬ 
dition, for very low supernova and magnetar energies, 
« IQSO ergs, the main light curve evolves 
slowly and a low luminosity double-peaked light curve 
can be seen. 

Three additional physical effects may further distin¬ 
guish the shock break out peak: 1) Though th e model 
light curves here and in Kasen & Bildsten (2010) assume 
a constant grey opacity, the true ejecta opacity likely in¬ 
creases inwards, given the higher temperature and ioniza¬ 
tion state of the inner regions, and the possible presence 
of synthesized iron group elements. A higher opacity in 



Fig. 7.— Approximate model light curves including both shock 
breakout emission (blue lines) and diffusive luminosity from mag¬ 
netar heating (red dashed lines). The magnetar parameters are 
Em = 5 X 10®^ ergs and tm = 2 days, corresponding to Bn = 
4, P = 2.5 ms. The supernova ejecta had kinetic energy Esn = 
1 X 10®^ ergs, mass 5 Mq, and opacity k = 0.1 cm^ g”'^. The right 
panel assumes 100% thermalization of magnetar energy; the left 
panel assumes inefficient thermalization according to Eq. |33| with 
Y = 0.1, A = 0.9. 



Fig. 8.— Approximate model light curves including both shock 
breakout emission and diffusive luminosity from magnetar heating. 
The magnetar parameters are Em = 5 X 10®^ ergs and tm = 2 days, 
corresponding to E 14 = 4, P = 2.5 ms. In the left panel, the 
supernova kinetic energy is Egn = 10 ®^ ergs and the mass is varied 
from Msn = 5 — 20 Mq . In the right panel, the supernova mass 
is Msn = IOMq and the kinetic energy is varied from Egn = 1 — 
3 X 10®^ ergs. The calculations assume 100% efficient magnetar 
thermalization and opacity k = 0.1 cm^ g~^. 

of the magnetar wind would delay the rise of the diffu¬ 
sive light curve; we discuss this point in more detail in 
the next section. 

4. MAGNETAR WIND THERMALIZATION 

If the magnetar wind does not thermalize efficiently at 
early times, the rise of the diffusive light curve will be de¬ 
layed, which will clarify the shock-breakout peak. There 
are physical reasons to think this delay in thermalization 
may occur in SLSNe. 

The spin-down luminosity of a magnetar is initially 
carried outwards by a Poynting flux. The magnetic field 
is initially strong close to magnetar surface, becoming 
dominated by its toroidal component outside of the light 


the interior would delay the diffusive light curve rela¬ 
tive to the shock breakout emission. 2) Deviations from 
spherical symmetry due to bipolar magnetar energy in¬ 
jection may lead to breakout occurring first along the 
poles, making the breakout emission more conspicuous 
from polar viewing angles. 3) Inefficient thermalization 


cylinder 

radius. Fventuall 

VI, reconnection inside the neb- 

ula (e.g. 

Forth et al. 

2013 

Mosta et al. 112014 

Bromberg 
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liovskoy|2(Jli; 

^ wil 

covert the magnetic energy to 

high energy particles. 

However, prior to the (uncer- 


tain) timescale for reconnection, the magnetar wind en¬ 
ergy is not inefficiently thermalized, and takes the form 
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of a magnetically dominated outflow that drives a shock 
through the ejecta. 

Even after reconnection dissipates the magnetic field 
energy, other physical effects may reduce the thermal- 
ization efficiency for a continued period of time. The 
dissipation of the wind energy at the termination shock 
or reconnection layers generates primarily high energy 
pairs. The injected pairs cool rapidly via synchrotron 
and inverse Compton radiation, producing high energy 
nieC^) photons that may in turn generate additional 
p airs by interacting wi th background thermal radia- 


where 


tion (Metzger et al. 2014). The optical depth to 7 — 7 
interactions is quantihed by the compactness parameter. 


TfnthO’T^’n 


2000L„,,45tday^^9 


(30) 


where T^m ,45 = ergs s“^, ug = Ht/lO® cm s“^, 

tday = f/l '^T IS the Thomson cross-section, Enth ~ 
Lnit is the approximate non-thermal energy injected into 
the nebula by the magnetar on the expansion time, « 
47 rr ^/3 is the volume of the nebula, and Vn ^ Vtt is the 
approximate nebula radius. 

If ^ ^ 1 then pairs produced by the first generation of 
photons upscatter additional seed photons to sufficient 
energies to create additional pairs. The details of this 
‘pair cascade’ are complex, but the net effect is to convert 
a sizable fract ion Y of the inj ected spin-down power into 
pairs (e.g., Svensson|l987). The total number of pairs 
N± in the nebula is set by the equilibrium between the 
rate of pair creation and annihilation 




YL^ 


rrieC^ 


and . 


(31) 


In equilibrium (iVj = 1V±), the Thomson optical depth 
of pairs across the nebula is 


■^es = = 


4ycrTT„ 


-I 1/2 




This equilibrium is reached on a timescale, teq — 
I6rn/3cr^g, which is short compared to the evolution 
timescale as long as S> 16i;t/3c. 

Most of the remaining fraction 1 — E of the energy 
released by the cooling pairs goes into a non-thermal 
power-law tail of radiation. The high scattering optical 
depth of the nebula traps these photons and delays 
their thermalization. On average, a hard photon must 
interact with the nebula walls ^ (I — A)~^ times before 
thermalizing, where the albedo A is the probability that 
the photon is scattered back into the nebula instead of 
being absorbed by the walls. The ‘lifetime’ of a hard 
photon is therefore bife = 1^(1 — A)where ~ (t"s -|- 
I)rn/c is the photon diffusion time required for a single 
nebula crossing. 

If tiife exceeds the expansion time, then non-thermal 
photons lose energy to adiabatic expansion before their 
energy can be thermalized. This reduces the effective 
rate of thermal energy pro duction to a fraction of the 
magnetar spin-down power (Metzger fc Piro||20I4) 


Lth = 


Tn 


1 + (tlife/t) ’ 


(33) 


tlife _ 

t c(l — A) 


rl/2 l/2,-l/2 
0-6 ^ _ ^^m.45'^9 ^day 


and we have assumed ^ 1 . 

Depending on the characteristic values of Y and A, 
suppression of the magnetar power due to thermaliza¬ 
tion can be import ant. A typical v alue of the pair mul¬ 
tiplicity is F ^ O.I ( Svenssoii||1987 ), although its precise 
value will depend on the nature of the pair cascade and 
deserves further study. The albedo depends on the ion¬ 
ization parameter, which sets the ratio of scattering to 
absorption in the e.iecta wall. Photoionization calcula¬ 


tions by |Ross et al.| (1999) (their Fig. 2) show a rather 
high albedo A ~ l).9 across a range of photon energies 
^1 — 30 keV, although for high photon energies inelastic 
scattering results in a higher absorbed fraction]^ 

For values Y = 0A,A = 0.9, thermalization of the 
magnetar wind will be suppressed for ~ 1 week following 
the explosion. Inefficient thermalization will not affect 
the dynamics of the shock, but will reduce the early dif¬ 
fusive luminosity from direct magnetar heating. We in¬ 
cluded this effect in our light curve calculations by using 
the suppressed magnetar luminosity Eq. [^to determine 
the supernova light curve. The results, shown in the right 
panel of Figure [^demonstrate that inefficient thermaliza¬ 
tion may serve to better distinguish the double-peaked 
light curve shape. 

5. DISCUSSION AND CONCLUSIONS 

The predicted luminosity and timescale of magne¬ 
tar driven shock breakout are compellingly similar to 
those seen in the double-peaked light curves of some 
SLSNe. As a concrete example. Figure compares ap¬ 
proximate light curve models to the observed bol ometric 
light curve o f SN-LSQI4BDQ (as constructed by Nicholl|| 
et al. (2015)). The models have Mgn = 25 Mq, Egn = 
5 X I0°^~ei^s, = 1-5 days, and = 5 x 10®^ ergs, 

although we do not claim these values to be optimal or 
unique. The adopted magnetar energy is high, but well 
within the possible r ange recently found for rapidly ro¬ 
tating neutron stars (Metzger et al.||2015). 

To see a clear double peak in the models of Fig re¬ 
quired that we assumed inefficient magnetar heating at 
early times. We tried two ways of impl eme nting this: 
the first used the suppression function Eq. the second 
simply assumed that heating was completely inefficient 
until a time t = 15 day, and 100% efhcient thereafter. 
The latter results in a clearer separation of the two light 
curve peaks, in better agreement with the observations. 
Clearly the uncertain details of magnetar wind thermal¬ 
ization are important in setting the precise shape of the 
early time light curve. The failure of the model to fit 
the observations at t > 250 days may also be due to 
a decrease in thermalization efficiency at late times (for 
distinct physical reasons). 


® Also note that we have assumed that photons absorbed by the 
ejecta are immediately thermalized. However, in reality if their 
energy is instead deposited in the electrons in the hot outer ionized 
layer (where the Compton temperature is much higher), additional 
Compton down scattering may be required to diffuse this energy to 
optical/UV wavelengths, resulting in an effective value of A which 
is even higher. 
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Fig. 9.— Comparison of the observed bolometric light curve 
of SN-LSQ14BDQ (black circles) to model light curves with pa¬ 
rameters Msn = 25 Mq^Esu = 5 X 10^^ ergs, £ni = 1-5 days, and 
Em = 5 X 10^^ ergs. Model 1 (blue lines) assumes that the magne- 
tar wind energy is inefficiently thermalized as described by Equa¬ 
tion |33| Model 2 (black lines) simply assumes that the thermal- 
ization is completely inefficient for £ < 15 days, and 100% efficient 
thereafter. The later approach provides more contrast in the peaks. 


The reasonable model fits shown in Figure suggests 
that the breakout scenario holds some promise for ex¬ 
plaining double-peaked SLSNe. However one should not 
put too much weight on this comparison (or the inferred 
physical parameters) given the number of approxima¬ 
tions that have gone into our model light curves. In 
particular, we have made coarse assumptions regarding 
spherical symmetry, the efhciency of thermalization, and 
the treatment of the radiative transfer. Detailed radia¬ 
tion hydrodynamical calculations are needed to make a 
more meaningful comparisons to data. 

If, in contrast to the models of Figure!^ thermaliza¬ 
tion is efficient at early times, the shock breakout peak 
is typically not easily distinguished from the diffusive lu¬ 
minosity from central magnetar heating. To explain a 
double peaked light curve in this case requires a spe¬ 
cific sets of parameters - a supernova with a large ejecta 
mass (Msn ^ 15 Mq) and/or a high explosion energy 
{Esn ^ 3 X 10®^ ergs). For high ejecta masses, one might 
expect the collapsing stellar core to form a black hole 
rather than a neutron star. Most of our discussion of 
driven shock emission applies equally to an engine p ow- 
ered by black hole accretion (Dexter & Kasen 2013), in 
which case the energy input is set by the rate of fallback. 

Deviations from our assumption of spherical symmetry 
presumably affect the luminosity and timescale of break¬ 
out. The actual geometry is likely bipolar, as even a 
sub-equipartition toroidal magnetic field can act through 
hoop stress to confine the nebular pressure along the ro¬ 
tational axis. The resultin g anisotropic stress may drive 
a weak, wide-angle “jet” (Bucciantini et al. 2007|) and 
the shocked ejecta will take the form of a broad “co¬ 
coon” that enshrouds that jet (as has bee n discussed 


for normal gam ma-ray bursts (GRBs), e.g., Lazzati & 


Begelman 2005). In SLSNe, the shock breakout of this 


“cocoon” would emerge continuously over a timescale of 
several days or longer, as set by the engine duration and 


the large size of the remnant. Our spherical analysis here 
may still be used to roughly estimate the dynamics, but 
with the input magnetar power enhanced by a factor of 
47r/D where D is the solid angle of the “jet”. 

Jet-like collimation of the energy should presumably 
result in a brighter shock breakout, at least for some 
viewing angles. This may make double-peaked light 
curves conspicuous even when inefficient magnetar heat¬ 
ing is not invoked. The light curve will depend on ori¬ 
entation, with the breakout emission being more promi¬ 
nent for polar viewing angles, and less so for equatorial 
views. Because the jet cocoon is expected to be fairly 
broad, bipolar, and non-relativistic, at least some break¬ 
out emission is likely to be emitted in all directions. 

We have not considered how ene rgy is thermalized 
behind the magnetar driven shock. Katz et al. (2010| 
show that, in supernova shock breakout, the post-shock 
gas and radiation reach equilibrium for shock velocities 
v/c ^ 0.1. The velocities expected here are much lower. 
Vs — Wej ~ 0.01c. Once the shock emerges, however, and 
if it continues to be driven into low density, optically 
thin ejecta, equilibrium may no longer be reached. In 
this case, one could look for some fraction of the shock 
heating emerging as non-thermal x-ray or radio emis¬ 
sion. Figure shows that, for typical parameters, the 
shock heating rates are ~ 10'^^ — 10"^"^ ergs s“^ at day 10, 
dropping to 10"^^ — 10^^ ergs s“^ by days 50-100. 

The dynamical effect of magnetar energy injection has 
additional observational consequences. An abrupt but 
small increase in the Doppler shifts of line absorption 
features (by an amount « 1000 km s“^) may occur when 
the photosphere recedes into the region of non-monotonic 
ejecta velocities created by the shock (see Figure[^. Fol¬ 
lowing this, the photospheric velocity should coincide 
with the motion of the swept up shell. For long spin- 
down times, the shell may still be accelerating at the time 
spectra are taken, such that the photospheric velocity in¬ 
creases with time, counter to the behavior expected for 
free expansion. Such an effect may have been observe d 
in the helium lines of SN2005bf ( Tominaga et al.| 

a dou ble-peaked Type Ib superriova that |Maeda __ 

(2007) modeled with magnetar heating. Once the mag- 
netar has shut off, the shell and photospheric velocity 
should approach a constant value over time. These ex¬ 
pectations may have to be modified to account for asym¬ 
metries due to an anisotropic magnetar wind or hydro- 
dynamical instabilities in the shell. 

The magnetar model is but one explanation of double- 
peaked SLSN light curves. While the first peak in SN- 
LSQ14BDQ was too brief and bright to be explained 
by ®®Ni heating, possible alternative mechanisms include 
cooling emission fro m a hyper-energetic SN explosion 
( Nicholl et al. 2015) or interaction with a dense GSM 
(Moriya & Maeda 2012 Piro 2015). In either case the 
necessary ejecta kinetic energy is large, 4 — 50 x 10®^ ergs, 
depending on the assumed radiu s of the star or GSM shell 
(Nicholl et al. 2015 Piro 2015). The second light curve 
peak requires a distinct niechanism, either central engine 
heating or interaction with an additional GSM shell at 
larger radius. While such a multi-component scenario 
can not be ruled out, it is appealing that the magnetar 
model alone may be able to reproduce the double-peak 
without introducing additional model parameters, or re- 


la et al. 
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quiring extreme values of the existing ones. 

Determining the fraction of SLSNe with double-peaked 
light curves would help discriminate the mechanism re¬ 
sponsible. Models that explain both peaks with CSM in¬ 
teraction require two distinct CSM shells - one low mass, 
nearby shell and one higher mass, more distant shell. 
There is no obvious reason why pre-supernova mass loss 
would conspire to frequently produce such a configura¬ 
tion. In the magnetar model, on the other hand, some 
early time emission from shock breakout is a generic con¬ 
sequence of the central energy injection. 

Spectra taken at the time of the first peak would 
also be diagnostic. Magnetar driven shock breakout is 
expected to produce a blue {Tes ~ 20,000 K) quasi- 
blackbody spectrum that is mostly featureless due to 
the high temperature and ionization state. Any de¬ 
tectable line features would be of high velocity (> 
10,000 km s“^). In the CSM interaction models, in 
contrast, one might expect narrow line emission from a 
photo-ionized, slowly moving CSM shell, or perhaps nar¬ 
row line absorption if a second, cold CSM shell exists at 
larger radius. 

While the magnetar model has been most frequently 
invoked to explain SLSNe, a lower level of magnetar pow¬ 
ering may occur in less luminous supernovae, perhaps 
in some cases producing a double peaked light curve. 
For example, for supernovae where the explosion energy 
and magnetar energy are both « 10®° ergs, the predicted 
shock breakout luminosity is only « 10"^^ ergs s“^, but 
produces a noticeable early peak in the light curve. 

If double-peaked SLSNe light curves are indeed due 
to shock breakout, this may indicate an interesting con¬ 
nection with the recently discovered class of very long 
duration GRBs. Greiner et al. ( 2015[ ) present observa¬ 
tions of a ^ 10^ s long GRB which had an associated 


super-luminous optical transient, both of which they ar- 

gue are powered by a magnetar. 
suggested that the jet in this even1 

Metzger et al. 

(2015I 

was just barel 

y pow- 


erful enough to escape the stellar remnant. For events 
with longer magnetar spindown timescales > 10^ s, or 
higher ejecta masses, a relativistic GRB jet may fail to 
emerge, however the underlying engine may still be re¬ 
vealed by the breakout of the magnetar driven shock (or 
“cocoon”) producing a double-peak optical light curve. 

While the calculations in this paper have outlined 
the main features of magnetar driven shock breakout, 
radiation-hydrodynamical calculations that include 
realistic opacities and, ideally, magnetic fields and 
departures from spherical symmetry, are needed for 
detailed predictions. The observational constraints on 
the magnetar model are now many: the timescale and 
luminosity of the first peak, the shape and brightness 
of the second peak, the photospheric velocity evolution, 
and the luminosity, color, and decline rate of the 
late time tail emission. Simultaneous fitting of all of 
these observables within the limited model parameters 
{B, P, Msn, Esn) constitutes a non-trivial test of the 
paradigm. 
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APPENDIX 


w hich h as frequently been applied 


To calculate approximate light curves, we us e a one-zone formalism 
to model magnetar powered light curves (e.g., Kasen & Bildsten 2010| |lnserra et al.||2U13p. To model the supernova 
light curve powered by direct magnetar heating we consider the evolution of the internal energy, Eint, of the bulk of 
the ejecta 

dEiat dV , , 

dt ~ ^ dt^ ^ ^ 

where V is the volume, p is the pressure. In the diffusion approximation, the radiated luminosity is 

. c d{Ei,,t/V) (Ei^t/V) 


Lsn = 47rR 


3k/9 dR 




Assuming homologous expansion (i? = vt) gives 

Esn ~ L/intt/tgjj where tsn — 


SKMsn 

Attvc 


R 


1/2 


( 2 ) 


(3) 


The timescale tsn is similar to the diffusion time td (Eq. although here we include the energy input by the magnetar 
to calculate the velocity, v = [2{Esn + E^)/Msn]^^'^. 

The formal solution to the differential equation Eq. is 

Esn = [ 2Li„(t'Asn)e“^‘ 


0 


(4) 
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The treatment here is clearly approximate, as it neglects the formation and expansion of the shell structure. Previous 
calculations have shown, however, that the one-zone formalism well reproduces the light curves from more detailed 
radiation-hydrodynamical models (Inserra et al. 2013). 

To calculate the luminosity due to magnetar driven shock heating, we solve an independent one zone model, usin g th e 
same integral expression Eq. I^but with the heating rate replaced with the shock heating rate Cgh (from Eq. |20[). 
The time dependent Cgh was aetermined from our numerical integration of the shell dynamics (see Figure]^. We 
further replace tsn with tbo (Eq. 16), as tbo gives the appropriate timescale when the diffusion time from the shock 
heated region equals the elapsed time. In these calculations we use Eq. to include the radiative loss term in the 
shell evolution, i.e. Lgn ^ 0 in the energy equation (Eq. 12). The two one-zone model light curves were summed to 
give the composite supernova light curve. 
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